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ABSTRACT 

The Lidov-Kozai mechanism ()Kozail [19621 : lLidov|[T"962l ) allows a body to periodically 
exchange its eccentricity with inclination. It was first discussed in the framework of 
the quadrupolar secular restricted three-body problem, where the massless particle 
is the inner bo dy, and later extended to the quadrupolar secular nonrestrict ed three 
body problem (|Harringtonlll969HLidov fc Ziglin|[l97a iFerrer fc Osacaij|l993 ). In this 
paper, we propose a different point of view on the problem by looking first at the 
restricted problem where the massless parti cle is the outer body . In this situation, 
equilibria at high mutual inclination app ear (iPalacian et aLll2006), wh ich correspond 
to the population of stable particles that IVerrier fc Evans! ( 20081 . l2009f) find in stable, 
high inclination circumbinary orbits around one of the components of the quadruple 
star HD 98800. We provide a simple analytical framework using a vectorial formalism 
for these situations. We also look at the evolution of these high inclination equilibria 
in the non restricted case. 
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1 INTRODUCTION 



As it is known, the secular three-body pr oblem aft e r nod e 
reduction ha s two d egrees of freedom (e.g. iPoincar 
iMalige et al.l I2OO2')). However, due to what iLidov fc ZiglinI 
( 19761 ) called a happy coincidence, this problem is integrable 



when it is expanded up to order two in the ratio of semi- 
major axes, i.e. at the quadrupolar approximation. Indeed, 
the argument of perihelion of the outer body does not ex- 
plicitly appear in the quadrupolar expansion of the secular 
problem, thus giving one more integral of motion linked to 
the eccentricity of the outer body. 

The limiting case where th e inner body has no mass 
has been extensively s tudied l|Kozail Il962l : iLidovl 1 19621 : 
iKinoshita fc NakaillioOTi ). We will call this problem the in- 
ner restricted problem, while the converse case where the 
two inner bodies are massive and the outer body is massless 
will be called the outer restricted problem. In the inner re- 
stricted case, the conservation of the normal component of 
the angular momentum enables the inner particle to periodi- 
cally exchange its eccentricity with inclination (the so-called 
Lidov-Kozai mechanism). The inner restricted model is well 
suited when the inner body has a small mass with respect to 
the other two. However, when looking at higher mass ratios, 
for example in triple star systems, this is no longer justified. 
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Since the Hamiltonian of the quadrupolar problem of 
three masses is very similar to that of the inner restricted 
problem when it is written in elliptic variables, the study of 
the massive problem has mainly focused on the dynamics 
of the two inner bodies (jHarringtonl Il969l : iLidov fc ZiglinI 
1 19761 : iFerrer fc Osacail Il994h . These previous works com- 
pletely classified the different dynamical regimes and bifur- 
cations, using the equations of motion of the inner binary. 



There is however another limit-case to th e massive prob- 
lem, w hich is the outer restricted problem. iPalacian et al.l 
(|2006l ) have studied this case and discussed the existence 
and stability of equilibria in the non-averaged system using 
the framework of KAM theory. We give here a very sim- 
ple model of the outer restricted case which provides an 
alternate formulation of these previous results and which 
is directly usable in an astronomical context. We also fully 
describe the possible motions of the bodies and give an an- 
alytical expression of t heir frequencies. We use this model 
to explain the results of IVerrier fc Evans! (|2008! . !2009l ), who 
find populations of particles at very high inclinations around 
one of the components of the double-binary star HD 98800, 
which are stable even under the perturbation of the other 
component. We then look at the quadrupolar problem of 
three masses from the perspective of the outer restricted 
problem and show how the inner and outer restricted cases 
are related to the general case. Vectorial methods as devel- 
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oped bv lBoue fc Laskaij (|2006l . l2009l ): lTremaine et all (|2009l ) 
are extremely well suited for this approach. 



Ml 
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2 SECULAR OUTER RESTRICTED PROBLEM 

2.1 Derivation of the Hamiltonian 

We consider here the case of a massless particle orbit- 
ing a central binary object. We do not restrict ourselves 
with respect to inclinations or eccentricities. The compo- 
nents of the binary have masses mo and mi, the binary's 
total mass is TVfoi = mo -I- mi and its reduced mass is 
/3i — momi/(mo -I- mi). The two massive components 
have barycentric positions uo and ui. We also denote S = 
mo/Moi and ri = ui — uo, and r2 is the position of the 
outer particle relatively to the barycentre of the inner bi- 
nary. Using these notations, the particle has the following 
Hamiltonian: 



H 



-G 



mo 



r2 - Uo 



+ 



mi 



r2 - ui 



(2.1) 



where f 2 = r2 is the canonical momentum associated to the 
position of the massless particle, r2. Since uo — —(1 — S)ri 
and ui = (5ri, we can rewrite the Hamiltonian as: 



2 



G 



mo 



mi 



r2 + (1 — S)ri\ |r2 — 3ri\ 



(2.2) 



We now suppose that n <^ r2 and expand the Hamil- 
tonian to order 2 in r\/r2- 
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r2 
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(2.3) 



The first term is the Keplerian energy of the particle 
interacting with the binary, seen as a point mass Moi. It 
is equal to — GA1oi/2a2, where 02 is the semi major axis of 
the particle. 

Since we are interested in the secular behaviour of the 
particle, we average this quadrupolar Hamiltonian over the 
mean anomalies of the binary (Mi) and of the particle (Af2). 
In order to do this, we first introduce four unit vectors: 
(i, j, k) are bound to the orbit of the binary, remain constant, 
and will provide a natural reference frame; w is bound to the 
orbit of the particle and will vary. More precisely, i points 
in the direction of the perihelion of the binary, k is colinear 
to the angular momentum of the binary, and j = k A i; the 
last vector w is colinear to the angular momentum of the 
massless particle. 

We can then compute the following averaged quantities, 
where quantities indexed with 1 relate to the binary, quanti- 
ties with index 2 relate to the particle, and u is an arbitrary 
fixed v ector (see for instance the appendix of lBoue fc Laskaj 
iiooi)): 



((ri.r2)') 



Ml 



2 I A/2 



f ('■i-(k.r2)^) 

+ ^(4(i.r2)^-(j.r2)^) 



^(l-ei)3/2 



(r2.u)^ 



(w.u)2 



,2 2ai(l-ei)3/2 • 

The substitution of these expressions in (|2.3|) yields 



(2.5) 
(2.6) 

(2.7) 



Mi,M2 



GMo 
2a2 



GPial 



8ai(l-ei)3/2' 



e? - i j + (k.w)^ 



ei(4(i.w) - (j.w) ) 



(2.8) 



Since the particle has no mass, the only variable ele- 
ment of the binary is its mean anomaly Mi. After averag- 
ing over this angle, it is no longer present in the Hamilto- 
nian. After averaging over the mean anomaly of the par- 
ticle, its semi major axis 02 becomes constant. Moreover, 
w — sin 12 sin 572 i — sin 12 cos 5^2 j + cos 12 k, so the argu- 
ment of pericentre UJ2 of the particle does not appear in the 
averaged Hamiltonian. Hence, at the quadrupolar order, the 
conjugate momentum associated to UJ2, i.e. the norm of the 
angular momentum of the particle G2 ~ \/GA4oia2{l — e^), 
is constant. Therefore the eccentricity 62 of the particle is 
constant. This fact is a feature of the quadrupolar expan- 
sion, not a property of the restricted problem. As such it 
remains true when the outer body has a non -zero mass (see 
section[3l. This is the happy coincidence that lLidov fc ZiglinI 
(1976) noted. Finally, only one degree of freedom remains, 
related to the couple (12,^22). 

If we drop the constant terms in (|2.8p , and introduce the 
mean motion ni of the binary into the Hamiltonian {n\a\ = 
GTWoi), we get the following expressiorQ (see also eq. 10 in 
(|Palacian et aLllioOQ )): 



{H) = [(k.w)^ - e?(4(i.w)^ - (j.w)^)] , (2.9) 



where 



3 { a^\'"'^ Pi 1 
a — —Til 



4 \a2j Moi (1-62)^ 



(2.10) 



This Hamiltonian can be rewritten in a very compact 
form as: 



(H) = —- 'w.T.w , where: 



/ -4ef 
T = aG2 ef 
I 1 



(2.11) 



(2.12) 



^ We will from now write (H) for the averaged Hamiltonian, omit- 
ting the subscripts Mi , M2 ■ 
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We also give the expression of the Hamiltonian using 
the inchnation and the node of the particle: 



{H) = — [2 cos^ i2 — e\ siv? 42(3 — 5 cos 20,2 



(2.13) 



2.2 Equations of motion 

As discussed in l|Boue fc LaskarlboOfil ) , the equations of mo- 
tion for w are simply obtained by: 



w = T^rVw (H) Aw 

(j2 



(2.14) 



After computing the gradient, we find: 



w = —a [(k.w)(k Aw) — ei(4(i.w)(i Aw) — (j.w)(j A w))] . 

(2.15) 

If we note x — (i.w), y = (j.w), and z = (k.w), we get 
the following system for {x,y,z): 



X = a(l - ei)yz ; 
ij = — a(l + 4ei)2;x 
i = baeixy . 



(2.16) 
(2.17) 
(2.18) 



In these variables, the fact that w is a unit vector and 
the energy integral translate into the following equalities: 



2 , 2 , 2 

X + y + z 

2 2 / J 2 2\ 

z -ej_(4:X - y ) 



h = Cst . 



(2.19) 
(2.20) 



The system of three differential equations (|2.16P " (|2.18p 
has thus two independent first integrals and is as such in- 
tegrable. It is also straightforward from these two relations 
that 



4ei ^h^l 



(2.21) 



2.3 Motion of a massless body around a circular 
binary 

In the case of a circular binary, the Hamiltonian and the 
equations of motion greatly simplifj|f|: 



(H) 



aG2 



k.w) 



w = — a(k.w)(k A w) 



(2.22) 



(2.23) 



The scalar product (k.w) — cos 12 remains constant, 
and the nodes of the orbit of the particle simply precess 
around the angular momentum of the binary, with a con- 
stant precession rate: 

A . 3 f aiV^^ I3i cos 12 , , 

fl = -acosi2 ^ ~-ni I — -Tl—Ti 2^7 (2.24) 

4 \a2j Moi (1 - ei)2 



^ The next non-zero term of the Hamiltonian which is the fourth 
order in (01/02) plays an imp ortant part in the circular c ase as 
has been discussed in detail bv lPalacian fc Yaneuaj (I2OO6I ). 






Figure 1. Intersections of the energy surfaces and the unit angu- 
lar momentum sphere (a) and its projection in the {x, y) plane (b) 
for ei = 0.5. Intersections of the energy surfaces and the unit an- 
gular momentum sphere (c) and its projection in the (x, y) plane 
(d) for ei = 0.2. 



This precession is equivalent to the precession generated 
by the quadrupolar potential of a circular and homogeneous 
ring of mass /3i and of radius a\ following an idea w hich 
can be traced back to Gauss (see (|Touma et al.ll2009i ) and 
references therein). 



2.4 Motion of a massless body around an elliptic 
binary 

2.4.1 Qualitative overview 

When the binary is elliptic, the situation changes and can- 
not be explained any longer by the quadrupolar torque of a 
circular ring. If we substitute z"^ in (|2.2Up using (|2.19|l . we 
get: 

(l-|-4e?)j;^ + (1 -e?)y^ = 1 - /i > , (2.25) 
+ y'^ + z^ ^1 . (2.26) 

The intersections of the energy surfaces and the nor- 
malized angular momentum sphere of the particle can thus 
be seen as the intersections of elliptic cylinders with the unit 
sphere. For a given value of the energy h, the extremity of 
the unit angular momentum vector of the particle w will 
move on the intersection of the corresponding cylinder with 
the unit sphere. Figures [T]a and c show these intersections 
for different values of the energy as dotted lines drawn on 
the unit sphere, in two situations where the binary has an 
eccentricity of 0.5 and 0.2 respectively. The three axes cor- 
respond to the scalar products x, y and z that are defined 
in section [23] 

There are four visible kinds of trajectories: closed tra- 
jectories around the two poles of the sphere {x, y, z) = 
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(0,0, ±1), and closed trajectories around the points 
{x,y,z) = (±1,0,0). 

When the extremity of the angular momentum of the 
particle w follows a trajectory around the north pole, it 
means that it precesses around the angular momentum of 
the binary k with an inclination that is strictly inferior to 
90°: in this case, the orbital motion of the particle is pro- 
grade relatively to the orbital motion of the binary. 

When the extremity of the angular momentum of the 
particle w follows a trajectory around the south pole, it 
means that it precesses around the opposite of the angular 
momentum of the binary, — k, with an inclination that is 
strictly superior to 90°: in this case, the orbital motion of 
the particle is retrograde relatively to the orbital motion of 
the binary. 

When the extremity of the angular momentum of the 
particle w follows a trajectory around one of the two points 
{x, y, z) = (±1, 0, 0), it precesses around the direction of the 
perihelion of the binary or the opposite of this direction. In 
this case, the inclination oscillates around ±90°. 



2.4-2 Frequencies 

The frequencies of these motions can be found analytically. 
Indeed, using equation (|2.25|l . we see that x and y are on 
ellipses or arcs of ellipses bounded by the unit circle (figures 
[T]b and d show respectively the cases where ei = 0.5 and 
0.2). Thus, there is an angle (p such that: 



h 



y = 



l±4e^ 



1 



Using (|2.19p . we can then express as: 



h + iel 



5e? 



h 



(2.27) 
(2.28) 

(2.29) 



1 ± 4ef 1 ± 4ef 1 - ef ^ ' 

There are two opposite values of z for each (j), reflecting 
the symmetry of the system with respect to the orbital plane 
of the binary. If we use expression (|2.29|) in combination with 
equation (|2.18p . we obtain the following equation for 0: 



± arcsin(l/fc), changes its sign (which is accompanied by 
a change of sign in the z variable), and the angle cj) librates 
between —(j>o and +(j)o- Thus, the projection of w on the 
orbital plane of the binary is bounded by the unit circle 
to stay on an arc of ellipse (|2.25|l . In this case, w precesses 
around the direction of perihelion of the binary, so that both 
the inclination and the node of the particle librate around 
±90°. 

In both cases, the period of the motion can be calculated 
with a simple quadrature using equation (|2.30|l : 



T - 



16 Mc 
3ni /3i 



7/2 



^{l-ei){h + 4ei) 



(2.32) 



where K{k^) is the elliptic integral of the first kind defined 
by: 



Jo 
Jo 



7r/2 



if it'' < 1 
if fc^ > 1 



(2.33) 



The last case where — 1 (or h = e\) corresponds to 
the trajectories that separate the previous two types. They 
link the points {x,y,z) = (0,±1,0), and the associated pe- 
riod is infinite. In the projection on the {x, y) plane, these 
separatrices form the ellipse which is tangent to the unit 
circle. Since all trajectories that are inside this ellipse corre- 
spond to the precession of w around k, the width Axsep of 
the separating ellipse in the {x,y) plane gives an indication 
on the proportion of such trajectories. Using equation (|2.27|) 
and the fact that h = e\ on the separatrix, we get: 



A:^scp — 2 



l + 4ef 



(2.34) 



Therefore, when the inner binary is circular, this width 
is equal to 2, the full width of the unit circle, and the only 
possible motion is precession of w around ±k. When the 
eccentricity of the binary increases, the width of the sepa- 
ratrix decreases to zero, which is a limit case since it can 
only be reached for a value of the binary's eccentricity equal 
to 1. The precession motions of w around ±i thus become 
predominant when the eccentricity of the binary grows. 



i = TQy(i^~e2)(^rM^yi"^ 



5e? 



1 - h 



1 - ef ft ± 4ef 



(2.30) 
(2.31) 



The constant fc^ is positive because of relation (I2.2ip . 
The value k = 1 defines a limit between two dynamical 
regimes. If fc'^ < 1, or equivalently if /i > ef, </> never van- 
ishes and the projection of w on the orbital plane of the 
binary moves along the full ellipse (|2.25|l . In this case, w 
precesses around the angular momentum of the binary, k. If 
z > Q the mutual inclination of the two orbits is always less 
than 90° so the orbital motion of the particle is prograde; 
conversely, if z < the mutual inclination of the two orbits 
is always superior to 90° so the orbital motion of the particle 
is retrograde. 

If A:^ > 1 (or h < e\), then vanishes for <jf>o = 



2.5 Comparison with numerical studies 

In IVerrier fc Evans! l2009l ). the authors investigate a fam- 
ily of particles at high inclinations around the binary HD 
98800 Ba-Bb, which remain stable even under the pertur- 
bation of an outer third stellar companion. They isolate a 
nodal precession imposed by the inner binary as the stabiliz- 
ing mechanism working against the destabilizing Kozai per- 
turbations of the outer companion. They run simulations of 
test particles orbiting the binary HD 98800 Ba-Bb using non 
secular equations. They observe the libration islands around 
12 = ±90° and 0.2 — ±90° that we discussed in the previous 
section. As they show their results in the (12 cos 0.2,12 sin Q.2) 
plane, we plotted the energy levels of the outer restricted 
Hamiltonian using these same coordinates for an easier com- 
parison. Figure[2]shows these levels for different values of the 
eccentricity. The c. panel in particular uses the same value 
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Figure 2. Energy levels of the Hamiltonian l|2.8| l in the 
(12 cos 02,42 sin O2 ) plane for values of the eccentricity of the bi- 
nary ei = (a), ei = 0.1 (b), ei = 0.79 (c), ei = 0.9 (d). 

for t he eccentricity of the bi nary (ei = 0.79) as figures 4 and 
5 of (jVerrier fc Evansll2009( ). 

Verrier and Evans notice no apparent structure in the 
dynamics of the (62,(^2) couple apart from the circulation 
of the perihelion. This is in agreement with the fact that 
the particle's eccentricity is constant at the quadrupolar ap- 
proximation. 

The authors also suggest that the projection of the an- 
gular momentum of test particles along the line of apses 
of the binary may be an integral of motion. From the re- 
sults of the previous section, it is straightforward to see that 
the projection x of the angular momentum of test particles 
along the line of apses of the binary is not constant. It varies 
with an amplitude that decreases to when the inclination 
of the particle approaches ±90°, which can be misleading 
when looking at numerical results for highly inclined par- 
ticles. However, the norm of the angular momentum of the 
test particles is an integral of the secular motion. 

Finally, the authors give a power-law fit of the period of 
the libration of the node with respect to three parameters: 
the eccentricity of the binary, the ratio of the semimajor axes 
02/01, and the mass ratio of the binary, 5. Their power- law 
is fitted using particles with fixed inclinations (85°). They 
give in their equation (5): 

/ X 3.37 

Tocer'-'^-^M^J . (2.35) 

By rewriting the mass dependences of equation (|2.32|) . 
we get the following analytical dependence with respect to 
the mass ratio and the semi-major axis of the binary: 

Tocf^j {5{l^5))-^ . (2.36) 

These two exponents compare very well with the fitted 
power law, in spite of the differences between the two mod- 
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Figure 3. Dependence of the period 1 12. 321 1 with respect to the 
eccentricity of the binary, in normalized units. The full line corre- 
sponds to the calculated perio d, while the dash e d line corresponds 
to the power-law fit given by I Verrier fc Evans! | |2009|) . We used a 
least squares method to fit the relative position of the two curves. 

els. The dependency with respect to ei is rather complex in 
equation (|2.32|) . and it is best compared in figure [3] 

The grid of initi al conditions for the particles in 
IVerr ier fc Evans"2009t) extends however from 3 to 10 AU 
for a binary separation of 1 AU, so the quadrupolar approx- 
imation may not be sufficient to fully describe the motion 
of the particles with the lowest semi major axes. In par- 
ticular, Verrier and Evans state that some low inclination 
particles show large eccentricity variations and even insta- 
bility. This could be due to a low initial semi-major axis and 
to resonances that are eliminated in our secular model by 
the averaging over the mean anomalies. 



3 PROBLEM OF THREE MASSIVE BODIES 

As we already stated, the quadrupolar secular three-body 
problem is still integrable when all the bodies have positive 
masses. As such, it is possible to show how the outer re- 
stricted problem we discussed in the previous section relates 
to the general c ase, and to t he in ner restricted case studied 
by lKozail (|l962l 'l and lLidovl l| 19621 '). We will first express the 
Hamiltonian of the secular quadrupolar problem using the 
same vectorial method as in the previous section in order to 
focus on the relative movements of th e orbits. In their s tudies 
of t he secular quadrupolar problem, iLidov fc ZiglinI (|l976h 
and iFerrer fc Osacan (|l994l ) have shown that this problem 
depends on two parameters. We will then point out which 
regions of parameter space are topologically equivalent to 
the outer restricted case, and which regions correspond to 
the inner restricted case, in order to show the continuity 
that exists between both situations. 

3 . 1 Hamiltonian 

Let us consider three masses mo , m\ and m2 , this time with 
m2 7^ 0. We note the bary centric coordinates and impulsions 
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(ui, Ui)i=o,i,2- As in the previous section, we suppose that 
the two bodies of indices and 1 form a binary and that the 
distance of the third body to this binary is much larger than 
the separation of the binary. We still note 5 = mo/(mo+mi). 
We first perform a canonical change of variables to Jacobi 
coordinates, 

ro = uo fo = uo + ui + U2 = (3.1) 

ri = ui " Uo fi = ui + (1 — 5)vl2 (3.2) 

r2 = U2 — 5uo — (1 — <5)ui f2 — U2 (3.3) 

I n these coord inates, the Hamiltonian of the three bod- 
ies is (|Laskarlll989l '): 



H 



2/3i 



Gm2 



2/32 
mo 



+ 



mi 



V2 — 5vi 



(3.4) 



\t2 + (1 - 5)vi\ 

where fi^^ = rriQ^ + mj""*", /3^"^ = (mo + m-i)^^ + m^"*", 
^1 — G(mo + mi) and ^2 ~ G{mo + mi + m2). 

Using the fact that ri <C r2, we expand the Hamiltonian 
to order two in ri/r2 as in the previous section: 



H 



A. 

2/3i 



Mi/3i 
ri 



+ 



A. 

2/32 



r2 



^ Pim2 /g (ri.r2)^ 



24 



(3.5) 



The first two terms are Keplerian energies and are equal 
respectively to — /ii/3i/2ai and — /i2/32/2a2, where ai and 02 
are the semi major axes of the inner and the outer body in 
our system of coordinates. 

We now average over the two mean anomalies Mi and 
M2 in order to get the secular part of the Hamiltonian. We 
will define 4 unit vectors which are analogous to the 4 vectors 
we used in the first section: (ii, ji, ki) are tied to the orbit of 
the inner binary; ii points in the direction of the perihelion 
of the inner binary, ki points in the direction of its angular 
momentum, and ji = ki A ii. The last vector k2 is colinear 
to the angular momentum of the outer body. In this section, 
the vectors tied to the orbit of the inner binary will no longer 
have fixed directions. 

Using the same averaging formulae as in the previ- 
ous section and using the fact that (ii.k2)^ -I- (ji.k2)^ -|- 
(ki.k2)^ = k2 = 1, we can write: 



{^)mi,M2 
1 



Mi/3i _ M2/32 _ 3_Gm2/3i 
' 2ai 2a2 8 (1 



s2)3/2 a: 



+ 2e\ + (1 - e?)(ki.k2)^ - 5e?(ii.k2)' 



(3.6) 



After averaging over the two mean anomalies, the semi- 
major axes are constant. There are thus four degrees of free- 
dom in the system, associated to the two eccentricities and 
the two inclinations. As we explained in the previous section, 
the argument of perihelion of the outer body does not appear 
in the quadrupolar expansion, and thus the norm of the an- 
gular momentum of the outer body G2 = P2 \//i2a2(l — 63), 
is constant. This implies that its eccentricity 62 is constant. 
Using the reduction of the nodes would leave only one de- 
gree of freedom in the reduced Hamiltonian, associated to 



the couple (ei,wi). The full reduction of the Hamiltonian 
and its expression in elliptical variables is the approach that 
has been used widely, since it yields a very simil ar Hamilto- 
nian function as in the inn e r restricted problem ( j Harringtonl 
1 19691 : iLidov fc Ziglinlll976l : iFerrer fc Osacailll993 i 

We want however to look at the motion of the nodes, 
or equivalently the motion of the vector k2 in the moving 
frame (ii,ji,ki) of the orbit of the second body. 

In order to easily compute the equations of motion, we 
introduce two vectors associated to the orbit of the binary 
that are colinear to ii and ki , a nd include in their nor m the 
eccentricity of the binary, as in (jTremaine et al.ll2009l ): 



1 - e^ki 



Ti = eiii 



(3.7) 



If we drop all the constant terms in equation 13.61 and 
use the above vectors, we get: 



iH)M„M. = - ^ [21? + (^i-k2)' - 5(Ti.k2)^] , (3.i 



where 



, 3 faiY^^ /3i 



and A^oi , ni are defined as in section [2] 
3.2 Equations of motion 

The components of ICi , Ii and k2 have the following Poisson 



irnpo 

bracketfl (Sorisov fc Mamaev.,2005. : Boue fc Laskar„2006 ): 



{lCii,JCij} = -^Y^K-ik , {k2i,k2j} = ~^^k2k , (3.10) 
Ai G2 

^ —^f^lCik , {lCii,Tij} = —-j^Tik , (3.11) 



Ai 



Ai 



where Ai — f3iy^fiiai and e^jfe is the Levi-Civita symboQ. 
The equations of motion for the three vectors are thus: 



Ai 



(/Ci A Vjci/Z-I-Xi A ViiJf) 



k2 = -7^k2 A Vk,^? . 

Lr2 



(3.12) 
(3.13) 
(3.14) 



In order to look at the motion of the vector k2 in the 
moving f rame (ii,ji,ki) of the orbit of the second body, 
we use as lBoue fc Laskail (|2006l ) the above system to derive 
equations for x = (k2.i1), y = (k2.ji), z — (k2.ki) and ei. 
Indeed, x = (k2.Ti)/|Ti|, z = (k2./Ci)/|/Cil, e-i — |Xi|, and 
y is obtained using the identity x'^ + y^ + z'^ = 1: 



We use the following convention, where pi are momenta and qi 
positions: {/,,} = ^ - ^ . 

^ ^ijfc = +1 if (itjtk) is an even permutation of (1,2,3), eijfc = 
— 1 is the permutation is odd, and e^j^ = in all other cases. 
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(3.15) 



y = — a'(l + Ae\)xz 
I G2 X 



5q fiisy + a 



= [(l-ei)(2-53;^)+5ei2^] (3.16) 



zxyz 



I G2 _ 
Ai 



(3.17) 



(3.18) 



The equations for x, y and z contain two terms: the first 
one is identical to the outer restricted system, and the second 
one includes the motion of the reference frame (ii,ji,ki) 
induced by the interaction with the third body. Note that 
when Ai is very large compared to G2 so that we can assume 
that G2/A1 is equal to zero, which corresponds to the case 
where m2 <C mo and mi, the above system is identical to 
the outer restricted system (|2.16|l - (|2.18p . 

The conservation of the total angular momentum C = 
G1 + G2 introduces the two main parameters of the problem. 
Indeed, 



A^ (1 - ) + G2 + 2Ai J 1 - elG2z = G 



(3.19) 



We note 7 = C/Ai, 72 = G'2/Ai. The above expression 
of the norm of the total angular momentum can be rewritten 
as a second degree equation giving -y/l — ef as a function of 
z using the two parameters 7 and 72: 



(1 - ei)' + 2722: Jl - ef + 72" - 7 







(3.20) 



The Hamiltonian can then be rewritten as: 
1 



{H) = --a'Ai72[z' + ei{2 - z' - 5a;')] . (3.21) 

The inequalities — 1 ^ z ^ 1 and ^ ei < 1 give the 
boundaries of the parameter space and the range of possible 
values of ei for any given couple of parameterqj (7,72): 



I7 - 72K 1 , 



I7 - 72I ^ e{ < min[7 + 72, 1] 



(3.22) 
(3.23) 



With these notations, the outer restricted problem of 
section [2] corresponds to the limit where 72 = 0, and in this 
case ei = ^1 — 7^ is constant as we saw. Note that when 
72 > 7, we have G2 > C, so this part of the parameter space 
contains only retrograde motions. Our aim in this paper is 
to show the continuity between the outer restricted case we 
studied in secti on [2] a nd the inne r restr icted case that was 
investigated bv iKozail (1962) and iLidovl (jT962). Both these 
problems lie in the region of parameter space where 7 > 72 
so we will restrict our study to this cas' 



^ The left part of the second inequality is strict if 7 = 72 . 
^ The other half of the parameter space (7 ^72) corresponds to 
retrograde motions which are of less physical interest and much 
more technical to study using our approach, in particular because 
equation 13.201 1 does not have a unique solution in this case. The 
interested reader will find a complete discu ssion of this case in 
llLidov fc ZiElinlll976l : iFerrer fc OsacaJll994l . 



In our case where 7 > 72, there is only one acceptable 
root to equation (|3.20l) . which is: 



-12Z + \l (722)2+72 - 72 



(3.24) 



This relation implies that ei is a growing function of z. 
Note that z = cos 42 , where 12 is the inclination of the outer 
body in the reference frame of the inner binary. As such, 
coplanar prograde motions (2 = 1) will always occur for the 
maximal value of the eccentricity of the inner binary: 



ei,, 



= Vl-(7-72)^ 



(3.25) 



Conversely, low eccentricities for the binary will be as- 
sociated to lower values of 2, and thus higher inclinations. 
Relation (|3.23p implies that the inner binary can only have 
a circular motion if 7 + 72 ^ 1. In this case, coplanar retro- 
grade motion (2 = — 1) is not allowed, and the lowest value 
of 2 is: 



20 



cos l2,max 



7 



7I 



1 



272 



(3.26) 



When 7 + 72 < 1 however, coplanar retrograde motion 
(2 = — 1) is possible and the associated value of the eccen- 
tricity of the binary is: 



ei, 



^1 - (7 + I2Y 



(3.27) 



3.3 Dynamical regimes 

In appendix [A] we briefly derive in the framework of the 
present study the fixed points of the system and the bound- 
aries of t he dynamical regimes in parameter space that are 
given m (jFerrer fc OsacaJl994l 'l. The fixed points are named 
as follow: the north pole is called A'^, and the south pole S; 
linearly stable fixed points are named E, as elliptic, and lin- 
early unstable points are named H , as hyperbolic; finally, 
signs are placed as indices to refer to the symmetry of the 
problem with respect to the two planes x = Q and y = G. 
There are three dynamical regimes in the region of param- 
eter space we study. 

In region O of figure UJ the parameter 72 = G2/A1 is 
small (less than 1/2). This is the case in particular when 
the mass ratio m2/m\ is small. Moreover, 72 -I- 37! < 1. The 
phase space is topologically equivalent to the outer restricted 
problem of section [21 The north pole, which corresponds 
to coplanar prograde motion with maximal eccentricity for 
the binary is linearly stable. There are two additional stable 
fixed points E± in the plane J/ = (see sect ion 1X4)1 . They be- 
long to the same family as the fixed points y — z = Q, x = ±1 
of the outer restricted problem that are re sponsible for the 
stable high inclination orbits observed by IVerrier fc Evanj 
(200& ). When 7 -I- 72 ^ 1, the south pole which corresponds 
to coplanar retrograde motion with minimal eccentricity for 
the binary, is also linearly stable. The a panels of figures 
and [6] provide a visualisation of the topology of this case 
When 7 -I- 72 > 1, the south pole is no longer accessible as 
stated in the previous section. It is however replaced by a 



^ Note that the south pole in figure [5] a corresponds to the 
out-most trajectory; this is an artifact of the coordinate map 
(12 cos £72i*2sinr22) which sends the south pole of the sphere onto 
the circle 12 = 180° . 
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Y 

Figure 4. Parameter Space. The dark gray areas are excluded 
by equation l|3.22|l . tlie liglit gray area correspond to tlie part of 
parameter space corresponding to 72 ^ 7 wliicli we do not study. 
The dotted line 7 + 72 = 1 separates the zone where there can be 
coplanar retrograde motion associated to a minimum eccentricity 
for the inner binary that is strictly higher than (below the dot- 
ted line) and the zone where the inner binary can be circular but 
the inclination is bounded (see section [A2| |. In zone O, the prob- 
lem is topologically equivalent to the outer restricted problem. In 
zones I and F it is topologically equivalent to the inner restricted 
problem, with zone I being equivalent to situations above the crit- 
ical inclination and zone F being equivalent to situations under 
the critical inclination. The letters a— f correspond to the values 
of the parameters used to plot the corresponding panels in figures 
[Hand 11 

stable trajectory at a maximal inclination given by equation 
(|3.26|) . as can be seen on panel b of figures [5] and |6l This 
trajectory corresponds to a circular inner binary (see sec- 
tion [^2|. Finally, there are two unstable points H± in the 
s = plane that belong to the same family as the unstable 
points of the outer restricted problem x = z = Q, y = ±1 
(see section . 

Panels c of figures [5] and |6] show the limiting case be- 
tween regions O and I. On this boundary, 7^ +872 = 1. The 
two unstable points H± are now located on the boundary of 
the accessible part of the sphere. 

Regions I and F of figure |4] are both in the part of the 
parameter space defined by 7^ + 872 > 1. In this zone, the 
problem becomes topologically equivalent to the inner re- 
stricted problem studied by Lidov (1962) and Kozai (1962). 
In the inner restricted case, there is a critical value of the 
inclination (cos 12 = -^/S/S) under which a circular inner bi- 
nary is always linearly stable, and above which a circular 
inner binary is always linearly unstable, giving rise to Kozai 
cycles. 

In region I of figure |4l the dynamical regime is topolog- 
ically equivalent to the inner restricted problem in the case 
where the inclination is superior to the critical value. The 
limit trajectory z = zo which corresponds to a circular inner 




-ISO' -90° Q' 90° ISO' -1S(1° -90' 0" 90° 1R()° 



Figure 5. Trajectories in the (12 cos ^2, i2 sinr22) plane for dif- 
ferent values of the parameters. See section 13.31 for a detailed 
discussion and appendix 1X1 for calculations, a; (7,72) = (0.8,0) 
outer restricted case with ei = 0.6; b: (7,72) = (0.8,0.25); c: 
7 = 0.8, 7I = (1/3)(1 -7^); d: (7,72) = (0.8,0.4); e: (7,72) = 
(1.08,0.4); f; (7,72) = (1.28,0.4). 



binary becomes linearly unstable. However, the north pole 
and the two fixed points E± are still stable. In the inner 
restricted phase space, when the inclination is superior to 
the critical value, there are two possible behaviours for the 
periastron of the inner particle: it can either circulate, or 
librate around ±90°. In our representation, the circulation 
case corresponds to trajectories around the north pole, and 
the libration islands correspond to the two fixed points E±. 
This is shown in panels d and e in figures [5] and [6] 

In region F of figure 3] the dynamical regime is topo- 
logically equivalent to the inner restricted problem in the 
case where the inclination is inferior to the critical value. 
Only one stable fixed point remains, on the north pole of 
the sphere, associated to prograde coplanar motion. This is 
shown in panel f in figures [S] and 

In both regions I and F, the parameter 72 = G2/A1 can 
take higher values. This is in particular true when the mass 
ratio 1x12/ m-i increases. 

The curve between regions I and F is linked to the crit- 
ical inclination that is defined in the inner restricted case. 
Indeed, along that curve, given by equation lAlOl we have 
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Figure 6. Trajectories on the unit angular momentum sphere for 
different values of the parameters. See section [3.31 for a detailed 
discussion and appendix lAl for calculations, a: {7,72) = (0.8,0) 
outer restricted case with e\ = 0.6; b: (7,72) = (0.8,0.25); c: 
7 = 0.8, 7I = (1/3)(1 -7^); d: (7,72) = (0.8,0.4); e: (7,72) = 
(1.08,0.4); f: (7,72) = (1.28,0.4). 



the following limits when 7 ^ 00: 



72 
7 



72 — 7 " 



(3.28) 



When G2 is very large compared to Gi, we can make 
the following first order expansion: 



72 — 7 = 



G2-C 



Ai 

G2- x/Gj + Gf + 2(G2.gT) 
Ai 

G2-G2(l + (G2.Gi)/Gi) 



Ai 



(G2.G1 

A1G2 



(3.29) 
(3.30) 
(3.31) 
(3.32) 
(3.33) 



As such, we see that along the border between regions 
I and r, when 7 and 72 both tend to infinity, we have the 
relation 



(3.34) 



Recall that z = cos 12 , where 12 is the inclination of the 
outer orbit in the reference frame of the inner orbit. Thus, 
the inclination of the inner orbit relatively to the outer orbit 
is ii = —22, and the above equation becomes: 



cos l-i 



(3.35) 



This relation is precisely the one giving the critical value 
of the normal component of the angular momentum of the 
inner body in the inner restricted problem. 



CONCLUSION 

We first studied the case of a massless particle orbiting a 
binary at a long distance, and, in the secular and quadrupo- 
lar approximations, gave a full analytical description of the 
motion along with the expression of the period of the sec- 
ular motion. When the inner binary is circular, only nodal 
precession takes place. However, when the binary is elliptic, 
libration islands appear at high inclinations, and these is- 
lands grow bigger when the eccentricity of the binary rises. 
IVerrier fc Eraing (|2008l . I2OO9I ) observe a similar nodal libra- 
tion in their study of the stability of particle populations 
in the quadruple stellar system HD 98800, and we showed 
that the analytical framework that we derived for the outer 
restricted problem is well suited to explain the results of 
Verrier and Evans. 

The quadrupolar secular three-body problem is still 
integrable when all the bodies have positive masses 
jHarringtonI Il969l : iLido v fc Zighn 1976; Ferrer & Osacar| 
1994). Using a vectorial for malism as ( Bouc & Laskar 20061. 
2009l : lTremaine et al.ll2009l ), we looked at this problem from 
the point of view of the outer restricted case. We showed 
how the outer restricted problem relates to th e gener a l case , 
and to the inner restricted case studied by iKozail (|l962l ) 
and lLidovll|l962h : when the mass of the outer body is small 
enough compared to the mass of the inner body, the gen- 
eral case behaves similarly to the outer restricted problem. 
When the mass of the outer body increases enough, the gen- 
eral case behaves like the inner restricted problem. We gave 
an expression of the boundary between these two regimes. 

The outer restricted problem and its generalization to 
the non restricted case provide an interesting starting point 
in the study of circumbinary planetary systems, such as the 
one discovered rec ently around th e eclipsing sdB-|-M sys- 
tem HW Virginis (|Lee et al.ll2009l ). In this system, the in- 
ner binary is very tight with a period of 2.8 hr, while the 
proposed planetary companions have periods of 9.1 yr and 
15.8 yr, so the quadrupolar expansion is fully justified. An- 
other field of application of the outer restricted problem 
is the study of the motion of stars o r biting around binary 
black holes (iMikkola fc MerrittI I2OO8I : iGillessen et~al] l2009l : 
iMerritt et aLll2009l '). 
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APPENDIX A: FIXED POINTS AND 
BIFURCATIONS 

The fixed points and boundari es presented in sectio n 13.31 
have already be en studied by iLidov fc ZiglinI (|l976l ) and 
iFerrer fc Osacarl ([l994f) . We briefly present here their deriva- 
tion in the framework of the present formalism. With the 
notations of section |31 we will limit ourselves to 7 > 72 . 



Al Poles of the sphere, x — y = 

This case correspond s to case 1 in section 5 of 
l|Ferrer fc Osacail 1 1994 ). Note that their sphere is con- 
structed using the eccentricity and perihelion of the inner 
binary, and is thus different from our angular momentum 
sphere. 

For all values of the parameters in the domain we study, 
the north pole z = 1, which corresponds to coplanar pro- 
grade motion, is a fixed point of the system. The associated 
eccentricity of the inner binary is: 



ei,, 



= Vl-(7-72)^ 



(Al) 



It is the maximal value of the eccentricity. This fixed point is 
always linearly stable. It is noted N in figures [5] and |6l Figure 
[5] shows the lines of equal energy in the (12 cos f22, 12 sin SI2) 
plane, and figure [§] shows these lines plotted on the sphere 
of unit angular momentum of the outer body = 1. 

When 7 -I- 72 < 1 (under the dotted line in figure |4]), 
the south pole z — —1 (noted S in the following figures), 
which corresponds to coplanar retrograde motion, is also a 
linearly stable fixed point of the system. The eccentricity of 
the inner binary is minimal and equal to: 



ei,: 



= n/1-(7 + 72)2 



(A2) 



Note that in this region of parameter space, the inner 
binary cannot be circular. 

When 7 -I- 72 ^ 1 (above the dotted line in figure |4]), 
the minimal eccentricity of the binary is as deduced from 
1)3. 23[) . The south pole z = —1 does not correspond to a 
real value of the eccentricity in this case. This limit how- 
ever is not a bifurcation strictly speaking. When crossing it, 
the stable south pole of the sphere is replaced by a stable 
trajectory at maximal inclination. 



A2 Circular Trajectories for the inner binary 

In the region of parameter space where circular trajectories 
exist for the binary (above the dotted line in figure |4]), the 
value of z which corresponds to such trajectories is minimal 
and equal to: 



zo 



2 2 

7 -72 ■ 



272 



(A3) 



The equations of motion on the small circle of the sphere 
z = zq are: 



X = Q'y(-2^o + 72(2 — Sx^)) , 

y = -a'x{zo + ■y2{2 - 5x^)) , 

z = , 

ex = . 



(A4) 
(A5) 
(A6) 
(A7) 



The right hand sides of equations 
for a certain value of x equal to: 



2 7' + 37I 
Xq = 



IO7I 



41 and (IA5I) vanish 



(A8) 



The curve 7^ -1-372 = 1 separates in figure|4]the regions 
noted O and I. We can distinguish three cases: 

(i) 7^ + 372 < 1. In region O of figure 31 < so there 
is no fixed point on the circle z = zq. Kb such, this circle is 
a trajectory for which the inner binary is circular and the 
outer orbit precesses at a fixed inclination given by ^2, max — 
ArcCoszo- Moreover, this trajectory is linearly stable. 

(ii) 7^^ -I- 372 = 1. There are two fixed points on the circle 
z = 20 at the coordinates (x = 0, y = i-^/l — Zq). 

(iii) 7^ + 372 > 1. In this case, we must also check that 
1 — a;o — Zq 0. The limit case where there is equality 

yields: 



vl 



572 - (4 107'')72' -f (57* - 87V 3) = 



(A9) 



This boundary limits the regions I and F in figure ID 
When solving the above equation for 7I and selecting only 
the relevant solution satisfying 7 > 72 , 7 + 72 ^ 1, we obtain 
a solution that corresp onds to equation 44 in section 5.2 of 
iFerrer fc Osacaij|l994 ) and that can be written using our 
notations as: 



2 

72 



2 -f 57^ - y6(V 



11 



7 + 72 > 1 



(AlO) 



In region I, yg > so there are four fixed points on the 
circle z = zq, at the coordinates (±a;o, ±J/o). They are noted 
H±±, in panels d and e of figures [S] and [S] Moreover, the 
trajectories that correspond to circular binaries are unstable 
in this zone. In region I' however, 1/0 < so we are again in a 
region of parameter space where there are no fixed points on 
the circle z = zo, and the trajectories associated to circular 
binaries are again stable. 



A3 The X ^0 plane 

When a; = 0, the only non trivial equation remaining in 
system (|3.15p - p.l8p is i = 0. Looking for a fixed point 
different from x — y = 0, we have to solve 



'l-e2 + 272 = , 
which after using relation (|3.20|l yields: 



ei 
x 



y^l - 72 - 37; 


± 



2 1 



z — — 



n/7^ + 37| ' 
272 

\/7^ + 37l ■ 



(All) 

(A12) 
(A13) 

(A14) 
(A15) 



We thus have two symmetric fixed points in the plane 
X = 0. They are noted H± in figures [5] and [6] For these 
fixed points to exist, the associated eccentricity must be a 
real number. As such, their domain of existence is the re- 
gion noted O in figure |4l This is case 2.1 in section 5 of 
jFerrer fc Osacaij|l994l ). 



High inclination orbits in the secular quadrupolar three-body problem 11 



These two fixed points are linearly unstable in their do- 
main of existence. Note that in the outer restricted problem 
(72 ~ 0), these fixed points are simply x = z — 0,y = ±1. 

A4 The y = plane 

When 3/ = 0, the only non trivial equation we must solve 
is 2) = 0. Here again, we look for another fixed point than 
X — y = 0, thus we have to solve: 

(l + 4e?)z+^2L=[(l-e?)(2-5s'')+5e?z'] = . (A16) 

Substituting 1 — in place of a;^ and then -^/l — Cj in 
place of z using (|3.20|) . we get: 

(1 - elf - (j' + + (1 - elf + ^(7^ - ,lf = 

(A17) 

This equation is t he same as equation number 40 in 
(I Ferrer fc Osacaij[l99i l. In our region of parameter space, 
there is at most one root which satisfies to the constraint 
l|3.23p . The curve separating the zone where there is one so- 
lution and the zone where there is no solution corresponds 
to the case where the limit value ei = is a solution, and co- 
incides with the boundary between regions I and V in figure 
|4] which is given by equation (|A10|I . 

When there is a solution, the value of ei can be trans- 
lated into a value of z using (|3.20p . Since y = Q, we get two 
values of a; = ±\/l — z'^, and there are thus two symmetric 
fixed points on the sphere, which are both linearly stable. 
They are noted E± in figures [5] and |6] When 72 = 0, these 
fixed points become simply y = z = Q, x = ±1, which are re- 
sponsible of the stable orbits at high inclination as discussed 
in the previous sections. 
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